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ABSTRACT 

We present a comprehensive analysis of the spatially resolved stellar population prop- 
erties of 26 early-type dwarf (dE) galaxies in the Virgo cluster. Using Lick/IDS absorp- 
tion line indices we derive simple stellar population(SSP)-equivalent age, metallicity 
and [a/Fe] abundance ratio. In particular, we focus on the comparison of the stellar 
populations between the central nucleus and the surrounding galactic main body. The 
stellar populations of the nuclei are, for most dEs, significantly younger than those of 
the respective galactic main bodies, with an average difference of 3.5 Gyr. We find only 
five dEs with significantly older nuclei than their galactic main bodies. Furthermore, 
we observe most dE nuclei to be more metal rich compared to their host galaxies. 
These age and metallicity behaviours are shown by almost all dEs brighter than Mr 
= -17 mag. 

The metallicity of both nuclei and galactic main bodies correlates with the total 
luminosity of the dEs. However, the metallicity of the nuclei covers a larger range 
(-1-0.18 to -1.22 dex) than that of the galactic main bodies, which all have sub-solar 
metallicity. The ages of dE nuclei show a statistically significant correlation with the 
local projected galaxy density within the cluster, such that younger ages are predom- 
inantly observed outside of the high-density central cluster region. The alpha-element 
abundance ratios are consistent with solar for both nuclei and galactic main bodies. 

We also examine the presence of radial gradients in the SSP parameters for a 
subset of 13 dEs (up to 1.2 kpc or 15 arcsec radius). We notice two different types of 
gradients, namely smooth profiles that include the nucleus, and profiles where a break 
occurs between the nucleus and the rest of the galaxy. Nevertheless, an overall trend 
of increasing age and decreasing metallicity with radius exists, consistent with earlier 
studies. The a-abundance ratio as function of radius is consistent with no gradient. 

Possible formation scenarios for the nuclei of dEs are discussed. The young and 
metal-enhanced population of nuclei suggests that these might have formed at later 
epochs, or the termination of star formation activity in the nuclei might have occured 
relatively late, perhaps due to continuous infall of gas into the central potential well. 
Our stellar population analysis suggests that the merging of globular clusters is not 
an appropriate scenario for the formation of most dE nuclei, at least not for the 
brighter dEs. We speculate that there might be different formation processes which 
are responsible for the formation of dEs and their nuclei depending on their luminosity. 

Key words: galaxies: dwarf - galaxies: evolution - galaxies: formation - galaxies: 
stellar content - galaxies: elliptical and lenticular, cD - galaxies: clusters: individual: 
Virgo 
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Early-type dwarf galaxies (dEs, Mb > -18) are the 
numerically dominant population in the present-day 
Universe (Sandage et al. 1985; Binggeli et al. 1987; 
Ferguson & Binggeli 1994). They also exhibit strong clus- 
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tering, being found predominantly in the close vicinity 
of giant galaxies, either as satellites of individual giants, 
or as members of galaxy clusters (Ferguson & Sandage 
1989). Although the dEs are characterized by their smooth 
appearance, having no recent or ongoing star formation 
and apparently no gas or dust content, the understanding 
of their origin and evolution remain major challenges for 
extragalactic astronomy. Stellar population studies show 
that dEs exhibit on average younger ages as compared to 
their giant counterparts, and also a lower metal content 
according to the correlation of metallicity and luminosity 
(Michielsen et al. 2008). However, past studies provided a 
wide range of ages (e.g., Poggianti et al. 2001; Rakos et al. 
2001; Caldwell et al. 2003; Geha et al. 2003; van Zee et al. 
2004), from as old as being primordial objects to dEs with 
recently formed young stellar populations. 

It appears that dEs themselves are not a homogeneous 
class of objects. Sub-structures such as stellar disks, faint 
spiral arms or bars are quite frequent among the brighter 
dEs (Lisker et al. 2006, 2007). Many dEs were found to 
contain a central surface brightness enhancement consistent 
with a point source on top of the galactic main body (e.g. 
Binggeli & Cameron 1991, 1993), referred to as so-called nu- 
cleated dEs. The studies from the HST/ACS Virgo cluster 
survey (Cote et al. 2006), with their high angular resolu- 
tion, not only verified the presence of such a distinct nu- 
cleus but also showed that nuclei are ubiquitous in bright 
dEs, covering a range in nucleus brightness. Interestingly, 
dEs with comparably faint nuclei that had not been iden- 
tified before Cote et al. (2006) show several systematically 
different properties as compared to dEs with bright nuclei 
(Lisker et al. 2007, 2008). 

Different studies of dE nuclei from difi^erent data 
sets found several contradictory properties for the nuclei 
(Grant et al. 2005; Lotz et al. 2004; Cote et al. 2006). Par- 
ticularly, the ground and space based data sets yielded 
different results. Grant et al. (2005) found that the nuclei 
are on average redder than their surrounding galactic main 
body. On the other hand, studies using HST observations 
(Cote et al. 2006; Lotz et al. 2004) measured the dE nu- 
clei to be slightly bluer than the galactic part. Furthermore, 
Cote et al. (2006), who used high quality data sets from the 
ACS Virgo Cluster Survey, proposed that the nuclei rather 
closely match the nuclear clusters of late type spiral galaxies 
in terms of size, luminosity and overall frequency. Another 
related scenario is also emerging: the recently discovered 
new (candidate) type of extremely small dwarf galaxies, the 
UCDs (Ultra Compact Dwarfs) with typical magnitudes of 
-13 <Mb< -11 (Hilker et al. 1999; Phillipps et al. 2001), 
might be the remnant nuclei of tidally stripped dwarf galax- 
ies (Bekki et al. 2003; Drinkwater et al. 2003; Goerdt et al. 
2008). 

The formation mechanisms of the nuclei of dEs are 
poorly understood and various possibilities have been pro- 
posed, also depending on the evolution and formation of dEs 
as a whole. As the nucleated dEs are preferentially rounder 
in shape, van den Bergh (1986) proposed that the nuclei of 
dEs could have formed from the gas that sank to the centre 
of the more slowly rotating objects. Since they predomi- 
nantly appear in highly dense environments, like the centre 
of a cluster of galaxies, the pressure from the surrounding 
inter-galactic medium may allow dwarf galaxies to retain 



their gas during star formation and produce multiple gener- 
ation of stars (Silk et al. 1987; Babul & Rees 1992), forming 
nuclei in the process. In both proposed scenarios the nuclei 
are formed along with the evolution of the galaxy itself, i.e., 
continuous star formation activity occurs at the dE cen- 
tre as time passes. Unlike that. Oh & Lin (2000) suggested 
that dE nuclei might have formed in a different way, namely 
through subsequent migration or orbital decay of several 
globular clusters towards the centre of their host dE. 

It is difficult to provide a definitive observational test of 
these different scenarios for nucleus formation. Nevertheless, 
we can gain some insight by comparing the different obser- 
vational properties, in particular relative ages and chemical 
enrichment characteristics, of the nuclei with their galactic 
main bodies, as well as with UCDs as their possible descen- 
dants. However, we need to bear mind that there may be a 
mixture of different formation scenarios. 

Our previous study based on this dataset (Paudel et al. 
2010, hereafter Paper I) has focused on the analysis of the 
inner stellar populations of dEs as a whole, without sepa- 
rating nuclei and galactic main bodies. Instead, our inten- 
tion was to see the variation of the inner stellar population 
properties with different morphological subclasses of dEs (cf. 
Lisker et al. 2007), using a much larger sample of Virgo dEs 
than in previous Lick index studies. We showed that dEs 
with different substructure properties (with/without disk 
features, Lisker et al. 2006) have significantly different stel- 
lar populations: dEs with disk features are younger and more 
metal rich than dEs without disks. Therefore we concluded 
that these dEs probably do not have the same origin, as 
they also differ in their distribution with local environmen- 
tal density in which they reside. By selection, all dEs in our 
sample contain a central nucleus, therefore it seems impor- 
tant to see the nature of the stellar populations of the nuclei 
and the surrounding galactic main bodies separately. And 
since there are different possibilities for the processes that 
form nuclei and also dEs themselves, we ask: can the nuclei 
thus tell us something about the formation history of dEs? 

This paper is organized as follows. In Section 2, we de- 
scribe the sample of Virgo cluster dEs, observation and data 
reduction in brief. In section 3, we describe the measurement 
of line-strength indices in the Lick/IDS system. Our main 
results from the stellar population parameters are given in 
Section 4 and are discussed in Section 5. Finally, we sum- 
marize our findings in Section 6. 



2 THE SAMPLE, OBSERVATION AND DATA 
REDUCTION 

Our sample comprises 26 nucleated dEs in the Virgo cluster. 
The sample properties such as position in the color magni- 
tude relation, total galactic luminosity, radial velocity and 
their local projected density within the Virgo cluster are 
described in detail in Paper I. The sample covers the full 
range of local density and includes the different morpholog- 
ical dE subtypes, i.e., 8 dEs with disks (dE(di)s) and 18 dEs 
without disks, which we hereafter simply refer to as dE(N)s. 
One dE(di) (VCC0308) contains a weak blue color excess in 
the centre, thus being referred to as a blue-centre dE (cf. 
Lisker et al. 2006). 

The observations were carried out at the ESO Very 
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Large Telescope (VLT) with the FORS2 instrument. The 
1" slit and 300V grism provide an instrumental resolution 
of ~11 A(Full Width at Half Maximum, FWHM). The other 
basic observational properties and the data reduction pro- 
cesses are described in detail in Paper I. 

We carefully checked the issue of scattered light during 
the reduction of the data, since the presence of a significant 
amount of scattered light could produce an artificial gradient 
in the measured line indices. Fortunately, our MOS-MXU 
setup utilized in this investigation provides the opportunity 
to quantify it. There are always free intra-slit regions where 
no light enters directly from the sky. After the bias sub- 
traction these regions should not contain any flux, unless 
scattered light were present. We thus calculate the average 
amount of light within such regions manually. We find that 
the mean is zero within the uncertainties, which are of the 
order of some hundredths of a count. The F0RS2 pipeline 
reduction produces the same result. It therefore confirms 
that there is no scattered light left in the spectra. 

In a different way, there is still the probability of mixing 
the nucleus light out to far beyond the central nucleus in 
case of bad seeing or instrumental blurring. To examine this 
effect, we also observed a star in an additional slit along 
with each target-field. Then, through the light profile of this 
star, we quantify the amount of such light at a radius of 
3" beyond the centre. Our measurements show that spread 
nuclear light is less than 1% of galactic light at 3" distance 
from the galaxy centre. The observed FWHM of the stars is 
always ~1.3" or less, consistent with this negligible fraction 
of starlight at 3" from the centre. 

2.1 Extraction of nuclear spectra and analysis of 
light profile 

Our goals in this paper are the measurement of simple stellar 
population (SSP) equivalent parameters (see Trager et al. 
2008) of dE nuclei and a comparison with the SSPs of the 
surrounding galactic main bodies. Additionally, if the signal- 
to-noise ratio (hereafter SNR) permits us, we wish to explore 
gradients in the SSP parameters, which helps to determine 
whether the SSP of the nuclei is very different from the 
rest of the galaxy or is just a continuity of a smooth SSP 
gradient at the centre of dEs. Although there is no pre- 
cise definition for what a nucleus is, the working definition 
used by several studies is that an excess of light from the 
smooth exponential (or higher order Sersic) profile of the 
rest of the galactic part is observed, looking like a compact 
source sitting at the centre of the galaxy. Because of its 
compactness, it is considered as a point source and repre- 
sented with a seeing convolved Gaussian light profile. Like- 
wise, the study of Grant et al. 2005 represents the nuclei as 
a point source convolved with Gaussian seeing. Cote et al. 
2006 used a slightly different approach, by fitting a two 
component core-Sersic model (Graham & Guzman 2003). In 
Fig. 1 we can clearly see for most dEs the change in the light 
profile at the centre (e.g. VCC0216, VCC0856, VCC0545, 
VCC1353 and VCC1945). On the other hand, VCC0308, 
VCC0990, VCC1261 and VCC1826 exhibit a rather smooth 
light profile. There may be several factors which produce 
such differences in the light profile even though all dEs in 
this sample are confirmed as nucleated from other photo- 
metric studies (Binggeli et al. 1985; Lisker et al. 2007). In- 



Table 1. Basic parameters and signal-to- noise ratio for our 
targets. 
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The second and third columns are the measured signal-to- 
noise ratio (SNR) per pixel at 5000A for the galactic-light- 
subtracted nuclear spectra and the galactic main body 
spectra, respectively. The fourth column gives the half-light 
semi-major axis in SDSS r from Lisker07. The fifth and sixth 
columns are total galactic and nucleus absolute magnitudes in 
SDSS r, applying a distance modulus of m — M = 31.09 mag 
(Mei et al. 2007), corresponding to d = 16.5 Mpc. The last 
column represents the amount (in fraction of total light of the 
central aperture) of light subtracted from the central nucleus 
spectrum. Nucleus magnitudes were derived as described in 
Paudel et al. (2010): a two-dimensional elliptical model image 
of the galaxy, based on a Sersic fit to the radial profile, was 
subtracted from the original image, taking into account the 
median SDSS PSF of 1.4" FWHM. The nucleus magnitude 
was then measured by circular aperture photometry with 
r = 2" centered on the nucleus; the error is estimated to be 
0.2 mag. 

sufficient spatial resolution or observed seeing which might 
blur the steeper light profile of the nuclei makes it harder 
to separate the galactic light profile. However, Cote et al. 
(2006) have observed the existence of a profile break in the 
case of VCC0856, VCC1261, VCC1355, VCC1407, VCC1661 
and VCC2019, reconfirming the existence of a nucleus at the 
centre of dEs with HST high resolution surface photometry. 

It is rather difficult to carry out an analysis of the stel- 
lar populations of nuclei alone, because the nuclei are al- 
ways situated on top of the underlying galactic main bod- 
ies. It is also hard to separate the galactic light from the 
central nucleus of such a faint object. The studies that have 
been done by Chilingarian (2009) and Koleva et al. (2009) 
provide results without galactic light subtraction from the 
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Figure 1. The light profile of dEs. The crosses represent the observed flux along the slit, and the solid line is the fitted exponential 
profile beyond 3" and extrapolated to the centre. 



nucleus. Although there is, in our sample, typically a fairly 
large domination of light from the nucleus as compared to 
galactic light at the photometric centre of the dEs, still a 
considerable amount of underlying light of the host galaxy 
can alter the observed properties of the nuclei. We therefore 
aim to reduce the galactic light contamination in the nucleus 
spectra, attempting a separate extraction of spectra for the 
nucleus and the galactic part. 

We extract the nucleus spectra from the central 0".75 
(i.e., 3 pixels). The region between 0".75 and 3" is not used 



for this extraction, to avoid any effects of nucleus light in the 
spectra of the galactic main body. We then integrate over 
the interval 3" to 8" from each side of the nucleus to ex- 
tract the spectra of the galactic main body (see Appendix A, 
Fig. Al). The individual spectra of galactic main body from 
the different side of nucleus were then co-added to produce 
a spectrum of higher SNR. In order to subtract the galaxy 
light from the nucleus, we determine a scaling factor by fit- 
ting the galactic main body's light profile (measured along 
the slit) by an exponential profile and extrapolating it to the 
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very centre, yielding the amount of galaxy light contained 
in the nucleus aperture (see Appendix A). Given the above 
considerations about the difficulty of separating nucleus and 
galaxy, we point out that our approach ensures the removal 
of a significant part, yet probably not 100% of galaxy light 
contamination. For those few cases where the central light 
profile looks rather smooth with ground-based data, and can 
only be disentangled with space-based photometry, our "nu- 
cleus" spectrum thus needs to be considered representative 
for the combination of nucleus and galactic central light. 

Before co-adding the spectra from the different sides of 
the galaxies, we analyze their slit profile to check for incon- 
sistencies or asymmetries, e.g. by contaminating objects on 
the slit. We find only one galaxy, VCC1945, has an asymmet- 
ric profile that deviates from a smooth exponential profile 
on one side. We noticed that a bright point source (fore- 
ground/ background or intra-galactic globular cluster) lies 
on one side of the slit. Therefore, we remove the spectrum 
from this side. For completeness, we also compare the spec- 
tra from the different sides of the galaxy before co-adding 
them, and we always find good agreement. Finally, the mea- 
sured SNR at 5OOOA for both the galaxy-subtracted nucleus 
spectra and the combined galaxy spectra is given in Table 
1. 



3 LINE STRENGTH MEASUREMENTS 

Before measuring the Lick absorption line indices from the 
flux calibrated spectra of the galactic main bodies and nu- 
clei, we also carefully checked whether any emission lines are 
present, particularly since some dEs show a fairly young nu- 
cleus. However, we do not detect any [OIII] emission, thus we 
do not correct the H/3 absorption for possible contamination 
by emission. If such emission were present, it would make 
the measured H/3 absorption smaller, and therefore derived 
ages older. On the other hand, it could be possible that we 
do not see any emission lines because of the low spectral res- 
olution. To quantify what strength of an emission line in a 
high-resolution spectrum (i.e, model of Vazdekis et al. 2010) 
would be smeared out in a low-resolution like ours, such that 
it is not recognized visually, we select a model spectrum of 
age 2 Gyr, and added an emission of H/3. We then degrade 
the spectrum to the low-resolution of 11 A. We find that 
the added emission line could have an effect of up to 12% 
on the measured absorption line strength, which reveals a 
relatively small effect on the age. 

Note that we have not applied a velocity dispersion cor- 
rection for the Lick indices, because the expected galactic 
velocity dispersion, (jgai ^ 50 kms~^, is significantly below 
our spectral resolution Oinstr ~ 280 kms^^. Therefore these 
corrections are not necessary. 

To measure the absorption line strengths from the 
spectra, we use the routine Indexf developed by N. Cardiel. 
It uses the definition of the Lick indices from Trager et al. 
(1998) and also derives the uncertainty in measured 
strength using Monte-Carlo simulations. Calibrations of our 
measured line strengths to the actual Lick system have been 
done as described in Paper I (Section 4.2 and Appendix B 

^ http:/ /www. ucm.es/info/Astrof/software/indexf/indexf.html 



in that paper). 

We use the method of Lick indices (Burstein et al. 1984; 
Worthey et al. 1994; Trager et al. 1998) as a tool for es- 
timating the stellar population characteristics. We trans- 
late our Lick index measurements into SSP-equivalent ages, 
metallicities, and a-element abundance ratios by compar- 
ing them to the stellar population models of Thomas et al. 
(2003) by "minimization, following Proctor & Sansom 
(2002). For this we use the nine indices RSf, H7F, Fe4383, 
H/3, Fe5015, Mg& , Fe5270, Fe5335 & Fe5406. Note that the 
SSP models assume all the stars were formed in a single 
burst and have the same age and metallicity. In fact, the 
galaxies may be a composite stellar system formed during 
several episodic star formation events, with different chemi- 
cal compositions in general. Therefore, our estimated stellar 
population parameters can be considered SSP— equivalent 
stellar populations. The correlation of age and metallicity in 
the model fitting is illustrated in Appendix B. 



4 RESULTS: AGES, METALLICITIES AND 
ALPHA-ABUNDANCE RATIOS 

In this section, we present the SSP-equivalent ages, metal- 
licities and a-abundance ratios of our sample dEs (Table 
2). Note that, in case of the least luminous dE, VCC0725, 
we find that the sky noise becomes dominant beyond the 
central aperture. Hence, we remove its galactic part from 
the sample and therefore provide no SSP parameters for the 
galactic main body of this dE. 

We can clearly see that the ages of the nuclei are sig- 
nificantly lower than the ages of the surrounding galactic 
main bodies (Fig. 2). The differences are more prominent in 
the disky dEs: only VCC1304 has a nucleus that is older 
than the galactic part. Moreover, we find that only four 
other non-disky dE(N)s (VCC1167, VCC1333, VCC1389 
and VGC1661 — see Sec. 2) have nuclei with significantly 
larger ages than the galactic main bodies. The median dif- 
ference in age between the galactic main bodies and nuclei 
is 3.5 Gyr. Examining Fig. 2 individually galaxy by galaxy, 
one can see that VCC0856 shows the largest difference (>10 
Gyr) in age between the nucleus and the galactic part. The 
nucleus of the blue centre dE VCC0308, while having a 
young age, does not show up as being special, having an age 
of 1.5 ±0.1 Gyr, similar to other dE nuclei such as VCC0216, 
VCC2019 and VCC1826. 

The metallicity distributions of the nuclei and the sur- 
rounding galactic main bodies also differ: the majority of 
the nuclei are relatively metal enhanced as compared to the 
galactic main bodies. However, it is remarkable that those 
nuclei that are older or equally old as the galactic part are 
also less metal rich than the latter. We find that the nu- 
cleus of VCC1308 has the highest metallicity of -1-0.16 ±0.12 
dex. For all dEs, the galactic main bodies have sub-solar 
metallicity. The a-abundance ratio from nuclei and galactic 
main bodies show a wide distribution. The nuclei of three 
dEs (i.e., VCC0308, VCC0389 and VCC0545) show signifi- 
cant a-enhancement as compared to their galactic part. On 
the contrary, four dEs, VCC0990, VCC0929, VCC1353 and 
VCC1861, exhibit a significantly enhanced a-abundance in 
the galactic part as compared to their nucleus. 
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Table 2. SSP-equivalent stellar population parameters for the nuclei and the galactic main bodies. 
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'^without subtraction of galactic light and does not have a measurement of SSPs from the galactic main 
body (see text). 



In the right part of Figure 2, the green vertical bars 
present, for comparison, the derived stellar population pa- 
rameters of the UCD sample of Evstigneeva et al. (2007). 
Note that we only use the published four indices (H/3, Mgb, 
Fe5270 and Fe5335). However, we use the same method of 
estimation for the stellar population parameters. The UCD 
ages and metallicities are consistent with old and metal poor 
stellar populations. Almost all UCDs have ages ~10 Gyr 
and metallicities vary between -1.25 to 0.13 dex. The [a/Fe]- 
abundances are always super solar in case of the UCDs, with 
a mean of 0.31 dex, which is 0.34 dex higher than the mean 
[a/Fe] of the dE nuclei. 

The relation between the stellar population parameters 
and the local projected number density of galaxies in the 
cluster is plotted in Figure 3. The local projected density 
has been calculated from a circular projected area enclosing 
the 10th neighbor. It seems that there is a correlation be- 
tween the local projected density and the ages of the nuclei. 
The Spearman rank order test shows a weak correlation of 
the ages and metallicities of the dE nuclei with the local 
projected densities. The correlation coefficients are 0.5 and 
—0.4, and the probabilities of the null hypothesis that there 
is no correlation are 0.2% and 4% for the age and metallic- 
ity, respectively. Unlike this, a similar test shows that the 
SSPs of the galactic main body do not have any relation 
with local projected densities. 

The relations between the stellar population parame- 



ters and the total galactic luminosity are presented in Fig. 
4. At the top of each panel, we also provide the trend of 
the differences in the SSP parameters between the galactic 
main bodies and the nuclei. It is clearly recognized that al- 
most all dEs brighter than Mr = —17 mag have younger 
and more metal-rich nuclei than the galactic main bodies. 
On the other hand, there is a relatively large scatter in the 
low luminosity region, and we can see that some of the nu- 
clei are as old and metal poor as the galactic main bodies. 
However, the sign of the differences in age and metallicity 
between galactic main body and nucleus are completely op- 
posite at the fainter and brighter end of the plot. As there 
exists a well-known metallicity-luminosity relation in early 
type galaxies (Poggianti et al. 2001), our sample also follows 
this relation for both nuclei and galaxies, i.e., the metallicity 
decreases with decreasing total galactic luminosity. The de- 
rived [a/Fe] values are fairly consistent with a roughly solar 
value for both nuclei and galactic main bodies. 

In the right panels of Fig. 4, we provide the number 
distribution (in the histogram) of stellar population param- 
eters of the nuclei (in blue color) and the galactic main bod- 
ies (in green color). It seems that the ages of the galac- 
tic main bodies have a bimodal distribution, but the small 
number of data points in each bin and the fairly large er- 
rors in the age measurement increase the uncertainty; the 
bimodality thus remains a qualitative impression. The age 
distribution of the nuclei is highly dominated by nuclei of 
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Figure 2. A comparison of stellar population parameters. Tlie SSPs from the different parts of tlie dEs are represented with vertical 
bars of different color: blue for the nuclei and red for the galactic main bodies. The faint background colors indicate the dE subtype: 
blue for the nucleated dE with disk and blue centre, green for the nucleated dEs with disks, and red for the nucleated dEs without disk 
features. The Virgo UCDs are represented by the green vertical bars with gray background. For the UCDs, wc used published values of 
line strengths from Evstigneeva et al. (2007) to derive the stellar population parameters (see text). 




• Galactic main body 

• Nucleus 





10 20 30 40 50 60 70 

Local projected density, sq.deg.~-' 

Figure 3. The age, metallicity and [a/Fe] versus local projected 
density. Green color represents the galactic main body and blue 
indicates the nucleus. 



younger ages. The metallicity distribution however appears 
much broader in case of nuclei than galactic main body. The 
nucleus metallicity ranges from slightly super-solar (+0.18 
dcx) to strongly sub-solar values (-1.22 dex), and interest- 
ingly all dE galactic main bodies have sub-solar metallicity. 



4.1 Stellar population gradients 

Due to the low brightness of dEs, it is always challenging 
to get spectra from their outer part with sufficient SNR to 
study stellar population gradients. Some attempts have been 
made to derive the stellar population gradients in the differ- 
ent cluster dEs (Chilingarian 2009 for Virgo, Koleva et al. 
2009 for Fornax). These studies used different methods to 
obtain SSP parameters, namely through spectral fitting with 
SSP models. Chilingarian (2009) observed either flat or neg- 
ative radial gradients in metallicity in his sample. However, 
due to the relatively high uncertainty in the age estimation, 
he did not draw conclusions on the radial behavior of ages. 
The study of Koleva et al. (2009) reconfirmed the result of 
the existence of negative metallicity gradients and found ra- 
dial age gradients in the dEs, with older ages at larger radii. 

In Figure 5—8, we present the radial profiles of SSP- 
equivalent age, metallicity and abundance ratio, measured 
in bins along the major axis of the dEs. It is interesting 
that we can divide these trends of SSPs in two groups. The 
first group are those dEs which exhibit a smooth trend of 
increasing the age and decreasing metallicity with radius, 
beginning from the nucleus, such as VCC0308, VCC0490, 
VCC0929, VCC1261 and VCC2019. In contrast, the second 
group shows a break in the SSP profile when going from the 
nucleus to the surrounding galactic part, with the latter hav- 
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Figure 4. The derived ages (top), mctallicities (middle) and 
[Q/Fe]-abundanco (bottom), plotted against r— band absolute 
magnitude (left). The blue eolor represents the nuelei and green 
color indicates the galactic main body. On top of each panel, we 
also show the difference in the SSP parameters, i.e. galactic part 
— nucleus. In the right panel, we provide the number distribution 
of the parameters. 



ing a nearly flat gradient, like e.g. for VCC0216, VCC0856, 
VCC1304 and VCC1355. 

Three dEs, VCC2019, VCC1261 and VCC0308, show 
a significant gradient in age and metallicity, having a rela- 
tively young and metal enhanced nucleus. Likewise, the ages 
of VCC0389, VCC0490, VCC0990, VCC0929 and VCC1407 
also seem to correlate with the radius. Our derived ages for 
VCC0856 agree with the result of Chilingarian (2009) that 
this galaxy has a flat distribution of ages beyond the central 
nucleus. In addition to that, we can also see such a flatness 
in the age distribution of VCC1355. VCC1261 presents the 
largest gradient in metallicity starting from slightly super 
solar down to a sub-solar value of —0.75 dex. Although we 
do not see any strong trend of [a/Fe] with radius in most 
of the cases, VCC0216 and VCC2019 display the opposite 
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Figure 5. The radial age profiles of selected dEs (here we select 
those dEs which have sufficient SNR at the last radial bin, 11" to 
15"). 



trend of decreasing and increasing of [a/Fe] with radius, re- 
spectively. 

In Figure 7 we show the age and metallicity distribution 
of our dEs in the different radial bins. Note that there is not 
always the same number of dEs in each radial bin: due to 
insufficient SNR in the outer radii for some dEs, those were 
omitted from the respective bins. The first 1" bin contains 
25 dEs, and the second, third and fourth bin contains 24, 20 
and 16 dEs, respectively. Therefore, the y-axis represents the 
normalized fraction in percent. It is easily noticeable that the 
distributions change with radius: the inner bin is dominated 
by young ages and shows a broader metallicity distribution, 
and the fraction of old ages and low mctallicities increases as 
we go outward, with the metallicity distribution becoming 
narrower. 



5 DISCUSSION 

In this paper, we have characterized the stellar population 
parameters from the different parts of dEs: the nuclei and 
the surrounding galactic main bodies. Our primary motiva- 
tion for this is to improve our understanding of the physi- 
cal mechanisms responsible for the formation of dE nuclei 
and the subsequent evolution of dEs themselves. As we now 
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in Fis. 5. 



discuss, our study makes two important contributions in 
this context: (i) to much more firmly establish the SSP- 
equivalent stellar population parameters of dEs and their 
nuclei (ii) to cast new light on the spatially resolved stellar 
population characteristics of dEs. 

The surrounding galactic main body is represented by 
extracting its spectrum from a 5" radial interval beyond 
3" from the centre, avoiding any contamination with hght 
from the nucleus. We expect that, due to our method of 
subtraction of the underlying galactic light from the nucleus 
spectra (see Appendix A), we obtained comparatively clean 
spectra of the nuclei, with the derived stellar population 
properties from such spectra well representing the nucleus 
stellar population. Nevertheless, as outlined before, in the 
cases of weak nuclei there is still a chance that the remaining 
galactic light contributes significantly, such that the nuclear 
spectra represent the combination of nucleus and "central 
galaxy light". To test for a possible bias due to this effect, 
we select those dEs which have galactic light fraction (see 
Table 1) larger than 50% at the central aperture such as 
VCC0990, VCC1353, VCC1355, VCC1407 and VCC1826, 
but all these nuclei have ages less than 5 Gyr, and agree 
fairly well with the average age of the nuclei in total. 

Generally speaking, stellar population gradients can be 
used as a proxy for the study of the evolutionary history of 
early type galaxies, since different formation models predict 
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different gradients. In a nutshell, monolithic collapse mod- 
els (Arimoto & Yoshii 1987) predict slightly steeper gra- 
dients than the hierarchical merging model (White 1980). 
These predictions, however, mainly apply to normal early- 
type galaxies (Es). In case of early-type dwarfs, different for- 
mation scenarios might be relevant, such as morphological 
transformation, or simply a primordial origin (also see the 
discussion in Paper I). Nevertheless, the overall distribution 
of age and metallicity at the different radial bins suggest 
that it occurs more frequently that the inner parts of dEs 
are younger and more metal enhanced than their outer parts, 
which is consistent with previous studies (Chilingarian 2009; 
Koleva et al. 2009). We also see two distict behaviours of ra- 
dial SSP profiles; the presence of flat proflles may be due to 
a particular galaxy structure (i.e., a faint underlying disk) 
or may be an indication of a different origin. Among the 
dEs with smooth SSP gradients, VCC0308 only has a very 
weak blue centre (Lisker et al. 2006), so it may well be that 
other galaxies have just a bit weaker colour gradients and 
were thus not labeled "blue-centre dE" previously. On the 
other hand, VCC0216 and VCC0856 have a similarly young 
nucleus as VCC0308, but not an age gradient in the galaxy 
itself, which might lead to having no colour gradient. 

Another key result emerging from our study is a very 
clear picture of the differences between the stellar popu- 
lations of the nuclei and the galactic main bodies of the 
dEs. To our knowledge, no spectroscopic study has yet per- 
formed such a comparison with a similar sample size. Stud- 
ies based on color differences (Durrell 1997, Cote et al. 2006, 
and particularly Lotz et al. 2004) flnd slightly bluer nuclei. 
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Figure 8. The radial profile of the a-abundance ratios of the 
dEs, selected as in Fig. 5. 



It is, however, not straightforward to interpret these color 
differences in the sense of stellar population properties, as 
we know that a degeneracy in the age and metallicity ex- 
ists with color (see also Appendix B). In contrast to the 
explanation of Lotz et al. (2004) of having more metal rich 
populations in the surrounding galactic main bodies, we find 
a metal poorer and older population in the galactic part on 
average. In addition to this, as Cote et al. (2006) note, there 
exists a color-luminosity relation for the nuclei. We also find 
that the metallicity of dE nuclei correlates with the total 
luminosity of dEs. 

We have seen that there is almost no correlation be- 
tween the ages of the galactic main bodies and the luminos- 
ity of the dEs. This might, at first glance, imply that the 
reason for the apparent age dichotomy in Paper I, finding a 
clear correlation with luminosity for the central stellar pop- 
ulations of dEs, was due to the nucleus contribution to the 
central aperture light. However, Fig. A2 of the Appendix, 
which compares the SSPs resulting from the nucleus spec- 
tra before and after subtraction of the underlying galactic 
light, actually tells us that this conclusion is not true: if the 
very central stellar populations of the galaxies, whose pure 
hght cannot be seen due to the superposed nucleus, would 
be so much older than the nucleus itself, the difference be- 



fore/after subtraction would be quite significant, which is 
not found. Instead, the figure tells us that the very central 
part of the galaxy does also reach, in most cases, almost the 
young age of the nuclei. Thus, in many cases it is really the 
age gradient within the galaxy that makes the galactic part 
surrounding the nucleus appear significantly older than the 
nucleus itself in Fig. 4. 

5.1 Evolution of dEs and formation of nuclei 

As we mentioned in the introduction, many studies have dis- 
cussed the origin of the nuclei of dEs together with the evolu- 
tion of dEs themselves. It is challenging to provide definitive 
observational tests of these different scenarios. Moreover, we 
argue in Paper I that not all dEs are the same class of object. 
The dichotomy in the age distribution of the galactic main 
bodies also supports the idea that one type of dEs may have 
a primordial origin (Rakos & Schombert 2004), being rela- 
tively old and metal poor. These might have suffered either 
early infall into the cluster potential or formed together with 
the cluster itself. The common idea is that internal feedback 
might be responsible for the removal of gas, with the con- 
sequence that star formation activity ceases at such early 
epochs. 

On the other hand, dEs with a relatively young and 
metal enhanced galactic main body likely have a differ- 
ent origin. As they are also preferentially brighter and of- 
ten host disk-structure, they might have formed through 
the structural transformation of a late-type spiral into a 
spheroidal system, triggered by the popular scenario of 
strong tidal interactions with massive cluster galaxies. Sim- 
ulations have shown that late-type galaxies entering in a 
rich cluster can undergo a significant morphological transfor- 
mation into spheroidals by encounters with brighter galax- 
ies and with the cluster's tidal field (Moore et al. 1996; 
Mastropietro et al. 2005). This scenario is unlikely to pro- 
duce the observed radial SSP gradients: either metallic- 
ity gradients must have formed in the late-type galaxies 
and somehow preserved during morphological transforma- 
tion (see the discussion in Spolaor et al. 2010), or accretion 
of leftover gas towards the centre of the galaxy would have 
to be responsible for the creation of such gradients. How- 
ever, the flat [a/Fe] profile implies a similar star formation 
time scale everywhere in the dEs. 

As we discussed above, the fairly different types of dEs 
with and without disk structure might have a different ori- 
gin. It is therefore even more difficult to explain the ori- 
gin of the nuclei of these dEs with a single scenario. How- 
ever, from this and previous studies, it is becoming clear 
that the majority of dE nuclei are unlikely to have formed 
through the merging of globular clusters: Cote et al. (2006) 
already explained the difficulty of this scenario with the 
luminosity differences, and additionally we find that most 
nuclei are fairly young and metal rich, at least in case of 
the brighter dEs (Mr ^ -17.25 mag). There are still the 
nuclei of some fainter dEs (i.e., Mr > —17.25 mag) which 
have fairly old and metal-poor populations, more resembling 
the stellar population properties of globular clusters. They 
might have formed through a different process as the nuclei 
of brighter dEs. 

The younger and comparably metal-rich nuclei support 
the idea that the central stellar populations of dEs were gov- 
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erned by continuous infall and accretion of gas in the centre 
of the potential well, building the nuclei. The brighter dEs 
also host disk features (e.g. residual spiral arms/bars) and 
these dEs themselves might have been formed through the 
transformation of late- type spirals (Sc-Sd types). High res- 
olution HST imaging has shown that such late-type objects 
frequently contain a compact nuclear cluster (Boker et al. 
2002, 2004), and Cote et al. (2006) observed that such nu- 
clear clusters have similar sizes to dE nuclei. Stellar popula- 
tion studies have shown that the majority of nuclear clusters 
have ages of few tens of Myr (Seth et al. 2006; Walcher et al. 
2006) with episodic star formation activity. Following the 
simplest interpretation, it could be that the present day dE 
nuclei are simply the nuclear clusters of the transformed 
late-type galaxies, and their star formation activity faded 
with the morphological transformation of the host galaxies. 
However, this scenario again fails to explain the observed age 
difference between the nuclei and galactic main bodies, since 
late type disks are also considered to host star formation ac- 
tivity throughout the inner region and disk. Alternatively, 
the truncation of star formation in the disk due to inter- 
actions could be more efficient than in the nucleus, which 
eventually leads to the development of age/metallicity gra- 
dients in dEs and makes the central nucleus younger and 
metal richer than the galactic main body. In any case, more 
detailed numerical simulations are required to test these hy- 
potheses. 

We find that dE nuclei exhibit fairly different stellar 
populations than UCDs. Particularly, the relatively older 
population (larger than 8 Gyr) and slightly super-solar a- 
abundance of UCDs may seem to create an inconsistency in 
the idea of dE nuclei being the progenitors of UCDs. Nev- 
ertheless, the current sample of UCDs is limited, and the 
fairly large spread in the stellar population properties of dE 
nuclei may allow the possibility of UCD formation in the 
Virgo cluster by the stripping of such dEs whose nuclei have 
old and metal poor stellar populations (Paudel et al. 2010). 
Therefore, a larger sample of UCDs and perhaps a more rig- 
orous comparison of SSP properties than this work is needed 
before any strong conclusions can be drawn. 



have a fairly large range (-(-0.18 to -1.22 dex). All galactic 
main bodies of the dEs have sub-solar metallicity. 

• While we see two distinct behaviours of SSP profiles 
(with and without a break) the overall trend of increasing 
age and decreasing metallicity with the radius is consistent 
with earlier studies. The a-abundance as function of radius 
is consistent with no gradient. 

• These observed properties suggest that the merging of 
globular clusters might not be the appropriate scenario for 
the formation of nuclei in dEs, at least not for the brighter 
dEs. The younger and comparably metal-rich nuclei support 
the idea that the central stellar populations of dEs were 
governed by continuous infall/accretion of gas in the centre 
of the potential well, building the nuclei. 

• The heterogeneous nature of the stellar population 
characteristics of dEs hints at different formation scenar- 
ios of dEs, similar to the conclusion of our previous study 
(Paudel et al. 2010). Our results suggest that the old, faint 
and metal-poor dEs are more likely to have a primordial 
origin, while those with relatively young ages and a higher 
metallicity and luminosity may have formed through mor- 
phological transformation. 
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6 CONCLUSIONS 

We have investigated the stellar population properties of 
the central nucleus and the surrounding galactic main body 
for a sample of 26 dEs in the Virgo cluster and compared 
the SSP-equivalent stellar population parameters of the dE 
nuclei with the ones of a small sample of UCDs. In addition 
to this, we have derived the radial profiles for age, metallicity 
and [a/Fe] abundance for 13 dEs. Our main findings can be 
summarized as follows: 

• We find that for most of the dEs the nuclei are signifi- 
cantly younger (~3.5 Gyr) and more metal rich (~0.07 dex) 
as compared to the galactic main body of the gala:xies. Only 
five dEs have significantly older nuclei than their galactic 
main bodies, and dEs with old and metal poor nuclei are 
more likely to be distributed in the dense region of the clus- 
ter than the dEs with young and metal-enhanced nuclei. 

• The metallicity of dE nuclei correlates with the total lu- 
minosity of dEs, and the observed metallicities of the nuclei 
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APPENDIX A: SUBTRACTION OF GALACTIC 
LIGHT FROM THE NUCLEUS 

The Figure Al provides a schematic view of the subtraction 
of galactic light from the central nucleus. First, we average 
the galaxy frame in the wavelength direction between 4000 
A to 5500 A, assuming that there is not any severe change 
in light profile from exponential with the wavelength. The 
fitting of the galaxy light profile with an exponential has 
been done only considering the galaxy light beyond the 3" 
from the centre, because we assume that the light from the 
nucleus should not be spread out to these distances, as the 
mean FWHM is 1.25" for our observations. 

The scaling of the galactic light to match the centre of 
the galaxy has been done by extrapolation of the light profile 
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Figure Al. A schematic view of tlic fitting of tfic liglit profile for 
VCC0490 and tlie binning processes. Tlie cross symbol represents 
tlie distribution of the observed total light (i.e galaxy + nucleus) 
and solid line represents the exponentially fitted light profile of 
the galaxy. The dashed line is the residual nucleus after the sub- 
traction of galaxy light, which represents the pure nuclear light 
profile. 
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Figure A2. The comparison of the SSP-equivalent parameters 
after and before subtraction of galaxies' light from the nuclei spec- 
tra. 



to the centre of dEs. The scale factor C has been calculated 
using the following equation, 



C 



i=-l 
~32 

E 

i = 13 



(Al) 



where is the flux from the best fitted galaxy profile (solid 
line in Fig. Al), and i is in pixel scale (i.e., 0.25") with the 
origin at the central peak of the observed slit profile of the 
galaxies. Then we subtract the galaxy light from the nucleus 



usmg 



(A2) 



here, F° is the observed hght in the frame. 

Although our exponential profiles of the galaxies 
are in good agreement with the observed profiles (see 
Fig. 1), some dEs have steeper profiles than exponen- 
tial (Janz & Lisker 2008) - VCC0389, VCC0929, VCC1167, 
VCC1254, VCC1348 and VCC1861 have n ?a2. Note, how- 
ever, this finding is based on fitting a much larger radial 
interval from the imaging data. In these cases, we again de- 
rived the galactic light profile for n = 2, which produced a 
better match for VCC0929. However, the calculated differ- 
ence of the amount of galaxy light which might be left at 
the centre when using n = 1 was less than 30% of the total 
central hght when compared to n = 2. Therefore, we always 
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Figure A3. The comparison of the SSP-equivalent parameters 
of the galactic main bodies (red), nuclei of dEs (blue), and the 
result for the combined central light from Paper I (black), where 
a central spectrum was analysed without separating nucleus and 
galactic main body. 



used the exponential profile for scaling the galactic hght to 
the centre for all dEs. 
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Table Al. Measured line strength indices from the nuclei of dEs after subtraction of galactic light and corrected to the Lick system. 
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"without subtraction of galactic light. 



Table A2. Measured line strength indices from the galactic main body of dEs (i.e., 3 to 8 arcscc radial interval) and corrected to the 
Lick system. 
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APPENDIX B: EXTRACTION OF SSP 
PARAMETERS 

It is well known that the age-metallicity degeneracy is a dif- 
ficult problem to estimate galaxy age and metallicity. How- 
ever, there are several different methods have been suggested 
to cope with this complication. By using the large number 
of indices and adopting the technique of Proctor & Sansom 
(2002), the effect of this degeneracy on the estimates of SSP 
parameters can be minimized. Fig. Bl, shows examples of 
the of Ax^ contours obtained with the method we have used 
to derive the SSP parameters, indicating the minimum with 
a diamond symbol. The contours are drawn with Ax'^ = 
2.3 (i.e., errors including 2 degrees of freedom (Press et al. 
1992, Section 15.6)). This shows that the typical la uncer- 
tainties we obtain on the SSP paramters are of the order 
of 0.1 dex. The effect of the age-metallicity degeneracy (e.g. 
Worthey et al. 1994) can be recognized in the tilt of the 
contours in the age Vs metallicity plot. 




0.2 

[a /Fe] , dex 

Figure Bl. Examples of Ax^ contours in different projection 
planes of age, metallicity and [a/Fe]-abundance space. 



